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ABSTRACT 

Thompson has argued that the Kozai mechanism is primarily responsible for 
driving white-dwarf binary mergers and so generating type la supernovae (SNe). 
If so, the gravitational wave signal from these systems will be characterized by 
isolated repeating pulses that are well approximated by parabolic encounters. 
I show that it is impossible to detect these with searches based on standard 
assumptions of circular binaries, nor could they be detected by analogs of the 
repeating-pulse searches that have been carried out at higher frequencies, even 
if these were modified to barycentric time as a function of putative sky position. 
Rather, new search algorithms are required that take account of the intrinsic 
3-body motion of the source as well as the motion of the Earth. If these eccentric 
binaries account for even a modest fraction of the observed SN rate, then there 
should be of order 1 pulse every 20 seconds coming from within 1 kpc, and there 
should be of order 200 detectable sources in this same volume. I outline methods 
of identifying these sources both to remove this very pernicious background to 
other signals, and to find candidate SN la progenitors, and I sketch practical 
methods to find optical counterparts to these sources and so measure their masses 
and distances. 

Subject headings: gravitational waves — white dwarfs — supernovae: general 



1. Introduction 



Except near final in-spiral, binary-star gravitational- wave (GW) sources are strictly 
periodic, and therefore a Fourier transform of the amplitude can be used to detect these 
sources. The only nuance is that one must also search over the sky by converting to a 
barycentric frame as a function of sky position. And, once a source is located, one can then 
"tune up" the signal by fitting for the phase and spatial orientation of the binary orbit. 

Until very recently, it was believed that the vast majority of individually identifiable 
binary sources would be on circular orbits. In this case the signal is sinusoidal, so a Fourier 
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transform is not mere ly convenie n t, it a lso contains a matched filter to the signal, and so is 



"optimal". However, iThompsonl (1201 if ) has shown that the great majority of white-dwarf 
(WD) binary sources probably have highly eccentric orbits. Within a "Fourier framework" , 
the resulting signal is represented by a discrete power spectrum, which covers the range 
[2n/P, Q], where P is the orbital period, (OP/27r) = (1 — e)~ 3 / 2 and e is the eccentricity. 

However, for eccentric orbits, Fourier transform no longer provides anything like a 
matched filter. For isolated binaries, one could in principle simply extend the search by 
including eccentricity and time of periastron as additional parameters. Then, since such 
systems are still strictly periodic (or almost so), one could fold the, say, 5- year GW data 
stream by the period (modulated by the annual motion of Earth as a function of putative 
sky position) to optimize signal detection. 



Unfortunately, the eccentric WD binaries predicted by IThompsonl (1201 if ) are, by their 
very nature, not isolated. The high eccentricities are induced by the presence of a third body 
with a separation that is ~ 10-100 times larger than the semi-major axis of the eccentric 
binary. The lower limit of this range is required for hierarchical stability and the upper 
limit to make the Kozai mechanism effective. Thus, several extensions of existing search 
techniques will be required to detect these objects. Note that the low mass density of WDs 
implies low-frequency waves, which are only detectable from space. 

Some intuition into this problem is gained by noting that the signal is completely domi- 
nated by pulses near periastron, and that during these passages the eccentric orbit is hardly 
distinguishable from a parabolic one. Here I develop this approximation and apply it to the 
problem to extracting science from gravitational pulse signals. 



The Parabolic Pulse 



Parabolic orbits may be parametrized by 
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Here, r = (x, y) is the separation between the two masses, b is the distance at periastron, 
%b is a parameter, Vl = a/ GM/b 3 , M is the total mass, and t i s time. In the qua drupole 
approximation, the gravitational strain tensor hij is given by ( jMisner et al.l Il973l MTW, 
Eqs. 36.3, 36.47) 
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where /|j aceless is the traceless part of the moment of inertia tensor Z^, and D is the source 
distance. 

Differentiation yields 

h xx 3y 2 r — Qx 2 b — 2r 2 b h yy 12xb 2 — 2r 2 b h zz 2b h xy y 3 
~A = 3^3 ' ~A = 3^3 ' ~A = ~3? ~A = ^3' ^ 

where 

bD Ml \0.1R Q J Vkpc/ v ; 

is the normalization of the GW amplitude, Mi and M 2 are the two masses, and M = Mi+M 2 . 
These waveforms are shown in Figure [TJ As expected, they are strongly concentrated is a 
small interval of time ±f2 _1 near periastron, where 

The total energy of the pulse is most easily calculated by taking the limit e — > 1 in the 
standard formula for the mean luminosity of an eccentric binary (MTW, Eq. 36.16a,b) 



85vr 2 „ M 
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This energy may be directly compared to the potential energy at periastron 
Bp„u« « 9 J VM o y VA/ S 7 VO.lfl^' 



Hence, if the iThompsonl (120111 ) Kozai mechanism is responsible for even one Milky Way 
supernova (SN) la per 500 years, then there must be one pulse per (500yr)/Q ~ 0.2 seconds. 
And even supposing that only 1% of these are within 1 kpc of the Sun, there would still be 
one such "loud" pulse every 20 seconds. Hence, disentangling these pulses is important both 
from the standpoint of understanding the sources and removing them as a background. 

In addition, if this mechanism is indeed important for generating SNe la, then it will 
almost certainly create a large population of in-spiraling WDs that lack sufficient combined 
mass to explode but still contribute to the gravitational pulse cacophony. 



3. Pulse Search Strategy I: Isolated Eccentric Binaries 

I begin by analyzing the simplified case of isolated eccentric binaries. As just noted, 
the real sources are not isolated, but this permits me to explore most of the relevant physics 
before introducing this additional complication. 
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Then because the pulses are periodic, they would in principle turn up in a Fourier- 
transform type search] 1 ] However, the signal-to-noise ratio (SNR) of such a search would 
be degraded relative to a matched-filter search by roughly {PVl/2-k) 1 ^ 2 ~ (1 — e)~ 3//4 . For 
example, for e = 0.99, the degradation would be a factor ~ 30, meaning that roughly 30 
times smaller signals would be accessible to a matched-filter search than a Fourier-transform 
search. 

Of course, there is some cost associated with enhanced sensitivity. For each period P, 
one would need to consider PQ different independent phases, where I adopt Q^ 1 = 40 s as 
the typical value. And one would have to step through the periods at AP = P/(f2T mission ), 
leading to a total (for each sky position) of A^ try = T mission P max f2 2 = 3 x 10 10 independent 
trials, where I adopt T mission = 5 yr as the lifetime of the mission and P max = 3 days as the 
maximum period that will be searched. And for each such trial, one would have to try several, 
perhaps 10, independent pulse widths Efficient algorithms have alre ady been developed 



for carrying out such searches in studies of b oth periodic GW so urces (lAbadie et al.l 12010 



and references therein) and transiting planets (IKovacs et al.ll2002l ). However, there is also a 
cost associated with the added risk of false positives. If one assumes Gaussian noise, then 
the minimum SNR required to avoid false positives is ~ ^J2 lnN try . If, for example, the 
number of independent trials goes from 10 4 to 10 15 , then the minimum SNR goes roughly 
from 4.3 to 8.3, i.e., a factor of 2. This is very modest compared to the factor ~ 30 gain in 
sensitivity due to matched filters. 



4. A Rough Estimate of LISA Source Counts 

To proceed further, I conside r the specific example of the Laser Interferometer Space 



Antenna (LISA. iPrince et al.ll2009l ). a proposed space mission composed of 3 antennae sepa- 
rated by L = 0.033 AU in an Earth-like heliocentric orbit. Lisa sensitivity peaks for sinusoidal 
periods at P ~ 150 seconds at a threshold of h ~ 10 -23 , with a full-width half-maximum 
(FWHM) of about 1 decade. Hence, it is extremely well-matched to pairs of 0.7 M & WDs 
with periastra b ~ 0.1 R Q , which have pulse widths corresponding to "sinusoidal periods" 
of about 2n/tt ~ 260 seconds. That is, for WD binaries with M = 1.4 M , the FWHM of 
LISA sensitivity corresponds to 

b = 0.Q7±° o f 4 R Q (FWHM). (8) 



1 Such searches must be conducted separately for every independent direction on the sky by transforming 
to the barycentric time for that particular direction. I will discuss searching in many directions below, but 
for the moment will assume one direction has been chosen and that time has been corrected to barycentric. 
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However, when comparing forecasted LISA sensitivity to Equation (T4|), one must adjust by 
a factor ~ (1 — e) _3//4 , to account for the fact that the signal only builds up once per period 
P, not once per 27r/fl Even so, for e = 0.99, LISA is sensitive to strains just 30 times higher 
than its nominal sensitivity, i.e., /i~3x 10~ 22 . Comparing to Equation and noting 
from Figure [1] that the signal peaks at a few times A, it is clear that LISA will be sensitive 
to such binaries to several kpc. 

The period distribution of sources is governed by a complex i nterplay of pre- W D and 



post- WD Kozai, binary evolution, GW emission, and tidal effects ([Thompson! 1201 ll ) . which 
are beyond the scope of this Letter. However, one can make a simple estimate by assuming 
that all pulses are the same. Then energy loss is governed by dE/dt ocl/P, which implies 
that the cumulative number of systems is N(P) oc P 1 / 3 . Hence, the population is weakly 
dominated by high-period systems, while pulse generation is dominated by tight systems: 
(dN/d\nP)/P oc P~ 2/3 . This is important. It means that most of the "noise" is concentrated 
in the highest SNR objects, so easiest to detect. But it also means that the most interesting 
sources are in the long-period tail. If we assume a typical "injection eccentricity" e max = 0.99, 
then 2/3 of the systems will have e > 0.9. And if the ensemble of sources is responsible for 1 
SNe la per 50,000 years within 1 kpc, then there are of order 200 such binaries within 1 kpc. 
Obviously, this is a very crude estimate, but the point is that there are potentially many 
more sources than there would be to supply the same SN la rate from circular orbits. 



5. Signal-to-Noise Ratio Estimates 



To make a more precise estimate of the expected SNR, I integrate the square of the 
profiles shown in Figure [TJ averaging over all Euler orientations, and restricting the integral 
to ±3fi _1 . I assume that the Lisa sweet spot permits detection of source s with rms strain 
h = A = 1.0 x 10~ 23 in T = 1 yr of observation at 5 ex fjPrince et al.lbood . pp 2, 21). I then 
find a SNR for a single pulse of 
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where F(Q) is the functional form of LISA sensitivity ( iPrince et al.ll2009l ) normalized to unity 
at maximum sensitivity at = 24 s. Note that this formula basically scales oc MiM 2 /D. 
The integrated SNR from summing all the pulses observed during the mission is larger by 
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For example, for Vt 1 = 40 s, and e = 0.99, we have P = 2ttQ 1 (1 — e) 3 / 2 = 3 days, 
so SNR = 25 SNRp U i se . Thus, even requiring a 9 a detection (see below) would permit 
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detections throughout the Galaxy, provided that multiple pulses from the same source could 
be so identified and thus "co-added". Moreover, individual pulses from within D < lkpc 
would be detectable. They, and their weaker cousins from farther away would constitute 
an incredible data mine if they could be interpreted, but a vast cacophony of noise if they 
could not. Since the great majority of this "noise" is due to the small number of binaries 
that have already been driven to shorter periods, it is important to note that for e < 0.1, 
SNR > 140SNR pulse , implying that these sources are much easier to extract and "remove" . 

Integrating the square of the time derivative of the wave forms in Figure [T] yields the 
timing precision of each pulse, a t — 0.8 fi _1 SNR~^ lse . Hence, if all the iVp U i se = T mission /P 
pulses from a single source can be successfully aligned, the angular direction precision from 
fitting to the correct barycentric pulse-delay pattern is (in radians) 



That is, sources within D < 1 kpc could be located within a few arcmin. Position measure- 
ments from LISA orbital motion alone would suffer an exact degeneracy in ecliptic latitude. 
The directional information relative to the spatial orientation of LISA is about AU/L = 30 
times worse than the orbital information, but it does not suffer from this ecliptic degener- 
acy. It would therefore be sufficient to break the orbit-based-direction degeneracy for nearer 
sources, but not more distant ones. 



As mentioned in the Introduction, the Thompson-Kozai WD binaries are all embedded 
in triples with semimajor-axis ratios c^/cii ~ 10-100. Hence 0.3 < g^/AU < 10. To illustrate 
the problems posed by this, I consider the case Mi = M 2 = M 3 = 0.7 M , where M 3 is the 
third body. Then, the binary will orbit the center of mass with an amplitude <22/3, which 
means that, if not corrected, the pulse signal will drift by ~ 170 s(a 2 /AU) sini, where i is 
the inclination. This is many times larger than the pulse width, especially considering that 
for a 2 < 6 AU, the system will complete at least half an orbit in 5 years. Hence, if the orbit 
around the third body is not included, the pulses will not align, even approximately, and the 
signal cannot be recovered. 

Hence, for each sky position (and its corresponding barycentric correction), and each 
combination of (Q, P, t per i), where t per i is the time of pericenter of the inner binary, one must 
conduct at least a 3-parameter search, corresponding to the period, phase, and amplitude of 
a sinusoidal orbit about M 3 . Then the number of such trials would be ~ [(a-2,max/3)fi] 3 ~ 10 4 , 




(11) 



6. Pulse Search Strategy II: Triples 
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where 6AU. For a 2 > a2,max, a simpler one parameter uniform-acceleration model 

would be adequate. 

Thus, the full search would be over An (AU f2) 2 = 2 x 10 3 sky positions, -P max r m i SS ion^ 2 ~ 
3x 10 10 inner binary periods and phases, [(a2,max/3)fi] 3 ~ 10 4 outer binary trials, and perhaps 
10 different pulse widths. This is 10 19 independent trials, which requires a detection threshold 
of SNR > 9. The sheer volume of computations is forbidding. As already mentioned, there 
are efficient algorithms for doing the 3 x 10 10 period search calculations. The problem is 
that these must each be done for 2 x 10 8 different configurations. Here I will simply assume 
that Moore's Law will handle this problem, although there is probably room for algorithmic 
improvements as well. In particular, the shorter-period binaries that dominate the "noise" 
will be detectable from subintervals T <C T mission , which will permit elimination of the outer- 
binary trials, and even allow standard Fourier techniques in many cases. Removal of this 
dominant "noise" will be essential to finding the more numerous high-eccentricity binaries. 



What parameters can be extracted from such observations, and how can the remaining 
degeneracies be resolved? I will argue that full resolution requires identification of optical 
counterparts. The counterparts that are easiest to identify are 3- WD systems, and these 
are also the most likely to yield key spectroscopic data leading to complete resolution. I 
therefore focus first on these systems. There will of course be a huge number of systems 
without counterparts, which could be subjected to statistical analysis, but the analysis of 
that problem is beyond the scope of this Letter. 

First, I review the observables. The best-fit barycentric correction gives the direction 
on the sky, and the waveform gives the three Euler angles of the binary (up to discrete 
degeneracies due to the quadrupole nature of GWs). These will be of use further below but I 
ignore them for the moment. There are also (at least) three orbital parameters for third body. 
The remaining four parameters that can be measured are A, Q, P, and t pcr i. The period 
is a parameter of fundamental interest since it helps classify the systems according to their 
evolutionary state and rate of progression toward WD collisions. Together with Q, the period 
also gives the eccentricity, which is again of independent interest and also enables a more 
precise waveform calculation (although I expect that this will be a very minor consideration 
for high eccentricity orbits). The time of periastron will also be important further below, 
but will be ignored for the moment. This leaves two observables 



7. Extracting Science 



A = 2 
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that combine four physical parameters, Mi, M2, 6, and D. It is important to keep in mind 
Q, P, and t per i are all robustly determined from distinct features in the data, while A is 
somewhat degenerate with the Euler angles, particularly the inclination and longitude of 
nodes. However, I ignore this degeneracy for the moment and assume that A is also well 
determined. Then, since there are four parameters and two measurements, there remain two 
degeneracies to be broken for which two additional pieces of information are required. Real 
additional "information" can only be derived from counterparts, but at the outset, one can 
gain a rough idea of b and D simply by assuming, e.g., Mi = M 2 = 0.6 M . Of course, such 
an assumption would make it impossible to derive the masses and thus determine whether 
the system was a viable SN la progenitor, which is arguably the most interesting science 
potential of the sample. But it would give a rough estimate of the distance. 

To make any further progress would require a catalog of WD candidates over the es- 
timated distance range of the sample. Optical identification of WDs would be extremely 
difficult if the third star in the system were a main-sequence star. However, a substantial 
fraction of tertiaries are likely to be descendants of intermediate mass stars, and so them- 
selves WDs. The most efficient wa y to construct a cata log of WDs (or multiple WDs) is a 



reduced proper motion survey (e.g.. lSalim fc Gouldll2003l ). This will be a natural by-product 



of the Large Synoptic Survey Telescope (LSST) for about 3/4 of the sky. Since WDs are 
typically My < 15.5, LSST sho uld reach several h undred pc or more, depending on its ex- 



act performance. For example, lAbell et al.l ( 120091 ) expect that at r = 24 (so D ~ 400 pc 



allowing for extinction A r ~ 0.5) the "main" LSST survey (1/2 sky) will achieve a proper 
motion precision of 1 mas yr _1 , corresponding to a transverse velocity error of 2 km s _1 . This 
is about 5 times better than is required to reliably identify a WD on a reduced proper motion 
diagram. Hence, even in the less- well covered northern 1/4 of the sky, LSST reduced proper 
motion diagrams should be adequate to about 400 pc. 

Recall from Section at these distances SNR > 60 and hence the position is known 
to < 5'. Since the surface density of WDs to this distance is only ~ 25deg -2 , there would 
be only a few candidate objects consistent with the position determined from GWs, even 
allowing for a factor ~ 2 uncertainty in the reduced-proper-motion distance estimates (and 
the smaller error in A for these high SNR objects). 

A spectrum taken near periastron would, by itself, positively identify the system as the 
LISA counterpart because the velocity difference between the two components would be of 
order lOOOkms -1 . Even, if the system turned out to be single-lined, the observed (brighter) 
component would be the less massive, and therefore would be moving at extremely high 
velocity. If the system were double lined, the radial velocity curve would give the mass ratio 
Mi/M 2 , while the spectra themselves would yield individual masses. The ratio of these could 
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be checked against the radial- velocity value, while the sum would yield sini, which could be 
checked against the value determined from the GW pulse profile. A similar exercise could 
be applied to the third WD, yielding a comprehensive picture of the entire system. Even if 
the more massive (so fainter) WD were beyond the detection limit, the degeneracies could 
be partly resolved by obtaining a trigonometric parallax. 

These determinations can probably be made for a large fraction of WD triples out to 
400 pc, which plausibly number in the dozens. Moreover, since their individual pulses can 
be detected, the algorithms to identify these sources are much simpler than those outlined 
above. 

Finally, I note that WD binaries with main-sequence companions can also be positively 
identified from the correspondence between the pulse timing residuals and the radial velocity 
curve of the companion. While the WDs would not be directly observable, precise astrometric 
measurements combined with a spectroscopic mass estimate of the companion would yield 
both D and M\ + M2, which (with Eq. [12] ) permit a complete solution. 

I thank Todd Thompson, Rubab Khan, and Scott Gaudi for useful discussions. The 
manuscript benefited greatly from careful review by the referee. Work supported by NSF 
AST-0757888 
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Fig. 1. — Four non- vanishing components of the gravitational waveform due to a parabolic 
encounter between two masses, which is the limiting case of a highly elliptical orbit. The 
periastron is on the x ("parallel") axis and orbital plane is determined by this and the 
"perpendicular" axis. These and the "vertical" components sum to zero at all times, since 
the matrix is traceless. To find the waveform seen at Earth, one would have to rotate 
according to the Euler angles and then project out only the transverse components. The 
wave-amplitude A and frequency Q are normalized to units of M , 0.1 R & and kpc, for 
masses, periastra, and distance. 



